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ABSTRACT 

We present a mechanism for the crystahine sihcate production associated with the formation and 
subsequent destruction of massive fragments in young protostehar disks. The fragments form in the 
embedded phase of star formation via disk fragmentation at radial distances > 50 — 100 AU and anneal 
small amorphous grains in their interior when the gas temperature exceeds the crystallization threshold 
of ~ 800 K. We demonstrate that fragments that form in the early embedded phase can be destroyed 
before they either form sohd cores or vaporize dust grains, thus releasing the processed crystalline dust 
into various radial distances from sub-AU to hundred- AU scales. Two possible mechanisms for the 
destruction of fragments are the tidal disruption and photoevaporation as fragments migrate radially 
inward and approach the central star and also dispersal by tidal torques exerted by spiral arms. As a 
result, most of the crystalline dust concentrates to the disk inner regions and spiral arms, which are 
the likely sites of fragment destruction. 

Subject headings: circumstellar matter — hydrodynamics — stars: protostars — stars: formation — 
accretion disks 



1. INTRODUCTION 

The embedded phase of low-mass star formation, start- 
ing from the formation of a disk and ending with the 
clearing of a parent cloud core, has recently drawn a 
renewed interest due to the possibility of planet forma- 
tion in this phase. Indeed, circumstellar d isks are most 
massive in the embedd ed phase (Andre ws fc Williamsl 
120051 : IVorobvovl l2QQ9al : iGreaves fc RicQ 120101 ) and nu- 
merical hydrodynamics simulations indicate that, un- 
der certain conditions, embedded disks can fragment 
and form either protoplanetary embryos or proto- 
brown dwarfs, especially in the disk outer regions 
where stellar irradiation and viscous heating are mer- 
ciful (e.g. Durisen et al. 2007; Stamatellos & Whitworth 
[2009; Bolev 2009; Machida et al.l l2010l : lVorobvov fc Basul 
HoiOb; Cha fc Nayakshin 2010). 

While disk fragmentation in the embedded phase seems 
certainly feasible, the prospects for the survival of frag- 
ments and formation of gas giants or brown dwarfs on dis- 
tances of order tens or hundreds AU are gloomy. Grav- 
itational instability in the embedded phase is strong, 
fuelled with a continuing infall of gas from a parent 
cloud core, and resultant gravitational and tidal torques 
are rampant. As a consequence, part of the fragments 
are disper sed by the tidal torques exerted by the spi- 
ral arms (IVorobvov fc Basul l2010bl : Holey et al. 2010). 
Others are quickly driven into the disk inner regions, 
and probably onto the star, due to the loss of angu- 
lar momentum caused by the gravitational interaction 
with the trailing spiral arms (jVorobvov fc Basul I2006L 
l2010bl : ICha fc Navakshin 2010T ). Only those fragments 
that happen to form in the late embedded phase, when 
gravitational instability and associated torques are get- 
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ting weaker, may open a gap in the disk and mature into 
gas giants on wide orbits but the pr obability for such 
a "for tunate" outcome is rather low (jVorobyov fc Basul 
I2010b[ ). 

Ironic as it may seem, but the destruction of fragments 
(hereafter, planet/brown-dwarf embryos or simply em- 
bryos) may be as important for the disk and stellar evo- 
lution as their survival followed by massive giant planet 
or brown dwarf formation. Embryos that are driven into 
the disk inner regions may form solid cores in their inte- 
rior via dust sedimentation or dense and compact atomic 
hydrogen cores via dissociation of molecular hydrogen. If 
one or both of those processes take place before the em- 
bryo is destroyed on its approach to the central star, a 
planet can emerge on orbit of order several AU or less 
(jNavakshinll2010al : iBolev et al.ll2010f ). If both processes 
fail, then the release of gravitational energy from the ac- 
creted embryo can trigger a luminosity outburst similar 
i n magnitude those of FU Ori onis or EX Lupi family 
(jVorobvov fc Basul 120061 l2010b[ ). In addition, short but 
intense bursts of accretion onto the star caused by the 
consumption of embryos can lower the lithium abundance 
in young solar- type stars (Baraffe fc Chabrier 20T3), ex- 
plain the luminosity spread in H-R diagrams of 1-Myr-old 
star clusters ( Baraffe et al.. .2Q09i) and resolve the long- 
standing luminosity problem (|Dunham et al.ll2010f ). 

On the other hand, embryos are sites of acce ler- 
ated d ust growth a nd hi gh-temperature process ing (Boss 
119981: iBoss et al.l 120021 : iBolev et all I2OIOI : iNavakshin 
l2010al |bl) and their in situ destruction by tidal torques 
may release processed dust directly into the disk outer 
regions. 

Whichever mechanism of embryo destruction is real- 
ized in nature, the consequences of this process are far- 
stretching. In this Letter, we report another interest- 
ing by-product of the embryo dispersal — the enrichment 
of embedded circumstellar disks with crystalline silicates 
formed in the depths of massive embryos via thermal 
annealing of amorphous dust grains. We use numeri- 
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cal hydro dynaraics siraulations of the disk forraation and 
long-term evolution complemented with equations de- 
scribing the transformation of micron-sized amorphous 
dust grains into the silicate form. We show that part 
of the crystalline-silicate-bearing embryos are destroyed 
before they vaporize dust grains in their interiors or form 
solid cores via dust sedimentation, thus enriching the gas 
disk with crystalline silicates on radial distances from 
sub-AU to hundred- AU scales. 

2. SOURCES OF CRYSTALLINE SILICATES IN 
CIRCUMSTELLAR DISKS 

Circumstellar disks originate during the gravitational 
collapse of rotating prestellar cloud cores. The degree 
in crystallinity of silicates in prestellar cores is similar 
to that of the interstellar medium, which has an up- 
per limit of 0.2 zb 0.2% in th e direction of the Galac- 
tic center (|Kember et al.|[2QQ3 ). On the other hand, the 
crystalline-silicate fraction in the inner regions of circum- 
stellar disks around young stellar objects (YSOs) can 
vary in wide limits f rom essentially none to almost 100% 
(|Watson et al.ll20Q9[ ). implying in situ production of crys- 
talline silicates in at least some YSOs. Two mechanisms 
for dust crystallization have been put forward. 

1. Evaporation of the original, amorphous dust grains 
followed by re-condensation und er conditions o f 
high temperature and density (e.g. lGrossmanlll972f ) 

2. Thermal annealing of the amorphous grains at 
temperatures (800-1300) K, somewhat below their 
vapo rization point, via viscous heating (e.g. Gail 
[20011 shock wave heating ( Harker fc Desch 2002). 
or disk surface heating during EX-Lupi-like out- 
bursts ( lAbraham et al.ll2009[ ). 

Both mechanisms are thought to work mostly in the inner 
1 AU from the star (except for shock heating that may 
operate at somewhat larger radial distances r < 10 AU) 
and some means of outward radial transport of crys- 
tallized dust is necessary to account for the non-zero 
crystalline-silicate fraction in comets at distances of or- 
der tens of AU (e .g. Gaili l200il iBockelee-Morvan et al.l 
[2002l : lGiesla][200^ . 

3. MODEL DESCRIPTION AND INITIAL CONDITIONS 

Our numerical model is explained in detail in 
I Vorobyov fc Basul (|2010bl ). Here, we provide only a brief 
overview and describe modifications done to the numer- 
ical code. We compute the gravitational collapse of ro- 
tating cloud cores in the thin-disk approximation and 
start our simulations from the pre-stellar phase, advance 
through the disk formation and early evolution phase, 
and terminate with almost a complete clearing of the na- 
tal core. The following physical effects are taken into ac- 
count in our model: stellar irradiation, background irra- 
diation with temperature Tbg = 10 K, viscous and shock 
heating, radiative cooling from the disk surface, and disk 
self-gravity. Turbulent viscosity is parameterized using 
the usual a-prescription. A spatially and temporally con- 
stant Oi is set to 0.005. 

In the current version, we have implemented a sim- 
ple mechanism describing the transformation of amor- 
phous dust grains into the crystalline form following the 



guidelines described in iPullemond et al.1 (|2006f ) . We as- 
sume that dust is passively transported with gas and 
neglect the dust radial drift caused by the dust-gas 
drag force. The latter simplification is of little con- 
sequence for the dynamics of dust grains with sizes 
< 10 [im on tin iescales of interest for the present pa- 
per, < 0.5 Myr (|Takeuchi fc Linll2002f ). Therefore, our 
results are strictly applicable only to sufficiently small 
dust grains showing sharp features in the 5-35 /im band 
of the mid-infrared spectrum. In this simple scheme, two 
additional continuity equations for the surface density of 
amorphous and crystalline silicates, Sa.s. and Sc. so I'e- 
spectively, are 

^^+Vp-(Sa.s.t;p) = -5, (1) 
^|^+Vp-(Se.s.t;p)=+>S, (2) 

where = Vrr-\-V(f)<j) is the gas velocity in the disk plane 
and S = i^cr^a.s.- All dust grains in a collapsing cloud 
core are initially amorphous, they crystallize at a rate Ucr 
provided by the Arrhenius formula 

Ucris-^] =2xl0^^exp(^-^y (3) 

where Tg is the gas midplane temperature, 2 x 10^^ 
is the latt ice vibrational fr equency, and Tc=38100 K 
([Bockelee-Morvan et al.|[2QQ2t) . We assume that dust and 
gas are thermalized. The initial gas and dust radial sur- 
face density profiles and angular velocities are given by 
equations (12) and (13) in Vorobyov fc Basu ( 2010b) and 
the initial dust-to-gas ratio is 0.01. The initial core mass 
is Mcore = 0.78 Mq and the ratio of rotational to gravita- 
tional energy of the core is /3 = 8.8x10"^. For simplicity, 
we also assume that all dust is made of silicate grains. 

A solution of the time-dependent part of Equation [2] 
(excluding advection terms) demonstrates that it takes 
about 40 yr to completely transform amorphous dust 
into the crystalline form at Tg = 800 K. The corre- 
sponding timescale for Tg = 900 K is 0.1 yr, while for 
Tg =700 K the crystallization timescale is 20000 yr. Ob- 
viously, the latter value is much longer than the typ- 
ical orbital period and crystallization temperatures of 
Tcr = 800 — 900 K are most relevant to circumstellar 
disks. Laboratory experiments tend to yield even higher 
Tcr > 1000 K (Hallenbeck et al. 2000), possibly due to 
the same timescale arguments. 

If temperature exceeds the evaporation threshold for 
dust grains Tgvap ^ 1400 K, we set both Ec.s. and Ea.s. to 
zero in the corresponding computational cell. The exact 
value of Tevap is density dependent and is calculated from 
the opacity tables of B ell fc LinI (|1994f ). 

Equations ([T]) and ([2]) are solved simultaneously with 
equations of hydrodynamics describing the formation 
and evolution of circumstellar disks (eqs (l)-(3) in 
IVorobvov fc Basul I20 10b). We employ the method of fi- 
nite differences with a time-explicit, operator-split solu- 
tion procedure in polar coordinates (r, 0) on a numerical 
grid with 512 x 512 grid zones. The advection part of 
Equations ([T]) and ([2|) (i.e., with S set to zero) is solved 
using a piece- wise parabolic advection scheme. The time 
update of Sa.s. and Sc.s. due to the source/sink term S is 
done using an (unconditionally stable) backward differ- 
ence scheme supplemented with a subcycling procedure 
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to guarantee accuracy, wherein the change in both Ha.s. 
and Sc.s. over one global hydro dynamical timestep is lim- 
ited to 10%. If this condition is violated, local subcycling 
with smaller timesteps is invoked until the desired accu- 
racy is reached. 

The radial points in our numerical grid are logarith- 
mically spaced. The innermost grid point is located at 
the position of the sink cell r^c = 6 AU and the size of 
the first adjacent cell is 0.08 AU. This corresponds to the 
radial resolution of Ar=1.3 AU at 100 AU. We impose 
a free outflow boundary condition so that the matter is 
allowed to freely flow through the sink cell to form a 
central star plus some dynamically inactive inner region. 

4. RESULTS 

The top row in Figure [1] presents the gas surface den- 
sity in the inner several hundred AU at three distinct 
evolution times since the formation of the central star. 
The middle and bottom rows show the corresponding 
gas midplane temperature and crystalline silicate frac- 
tion ^c.s. = Mc.s./(Mc.s. + Ma.s.), respectively, where 
Mc.s and Ma.s. are crystalline and amorphous silicate 
masses. The red circle in the coordinate center repre- 
sents schematically the central star plus sink cell. The 
three time snapshots in Figure [1] are chosen specifically 
to illustrate disk fragmentation and embryo formation 
(red dots at 50-100 AU). We note that each snapshot 
shows in fact different embryos — their inward migration 
and/or tidal destruction timescales are of order several 
orbital periods (Torb = 700 yr for M* = 0.25 M©, 
r = 50 AU, t = 0.074 Myr and Torb = 1600 yr for 
M* = 0.4, r = 100 AU, t = 0.178 Myr), which are con- 
siderably shorter than the time span between the images 
(- 5 X 10"^ yr). 

We find that the gas temperature Tg in the interiors 
of massive embryos can exceed 800 K (see middle row in 
Figure [1]) , necessary for the thermal annealing of amor- 
phous dust grains to commence. The prospects for disk 
enrichment with crystalline silicates stored in these em- 
bryos will depend on characteristic timescales for inward 
migration and/or tidal destruction of embryos (rdest), 
dust sedimentation and solid core formation in the em- 
bryo interiors (rged), and dust vaporization at Tg > Tevap 
('Tevap)- If 'Tdest is Smaller than both r^ed and Tevap, then 
all processed dust will be released into the disk at various 
radial distances, depending on which destruction mecha- 
nism prevails. As current numerical simulations indicate, 
embryos may be destroyed by tidal torques from spiral 
arms, preferentially releasing the crystalline silicates to 
large radial distances where these embryos form. Al- 
ternatively, embryos may be photoevaporated and torn 
apart by stellar irradiation and tidal torques as they mi- 
grate radially inward and approach the central star. In 
this case, part of the processed dust will be brought into 
the disk inner few AU (through the sink cell) and prob- 
ably onto the star and part will be pushed to larger disk 
radii during a transient episode of disk expansion (caused 
by the conservation of angular momentum) following the 
embryo tidal disruption. It is thus feasible to deposit 
crystalline silicates to various orbital distances starting 
from sub-AU scales and up to several hundred AU, ac- 
counting for the crystalline silicate features observed in 
the spectra of the Solar system comets. This process of 
crystalline silicate enrichment is illustrated in the bot- 



tom row of Figure [T] — there is a notable increase in ^c.s. 
with time. Crystalline silicates concentrate in the disk 
inner regions and at /near spiral arms, which are likely 
sites of embryo formation and destruction. This implies 
that the abundance of crystalline silicates may increase 
with radius, at least in the early disk evolution, if most 
of the embryos are tidally destroyed in the disk outer re- 
gions rather than driven into the star. If such an increase 
is observationally confirmed, it could become a hallmark 
for the proposed process. 

On the other hand, if Tdest greater than either Tged or 
Tevap, then most of the processed crystalline dust will 
be either locked up in solid terrestrial-like cores or va- 
porized. The latter mechanism is implemented in the 
current numerical model, while the effect of the former 
is discussed in more detail later in the text. 

To better illustrate the process of disk enrichment 
with crystalline silicates. Figure [2] presents various model 
characteristics as a function of time. In particular, the 
crystalline silicate fraction ^c.s. and amorphous silicate 
fraction ^a.s. = Ma.s. /(Mc.s. + Ma.s.) are plotted in pan- 
els A and B, respectively. More specifically, black, green, 
and blue lines depict the corresponding fractions in the 
disk (^^i«^ and ^^If ), envelope (C^s^ and ^I'^l), and sink 
cell (^^^g and respectively. We distinguish between 

the disk and envelope using a scheme based on the typ- 
ical transitional gas surface densit y and radial g as ve- 
locity field explained in detail in .Vorobvovl (|2010f ). We 
note that ^^'^^ and are in fact crystalline and amor- 
phous silicate fractions stored in both the disk and the 
embryos. Panel C shows the crystalline silicate fraction 
in the disk (black line) and maximum midplane gas tem- 
perature T^^^ (red line), whereas panel D presents the 

mass accretion rate onto the star M (red line) and the 
crystalline silicate fraction in the disk (black line). 

There are seven episodes of crystalline silicate produc- 
tion in the disk manifested by a sharp increase in ^^s^^' 
all associated with the formation of embryos massive and 
hot enough to crystallize dust in their interior. Episodes 
4, 6, and 7 (from left to right) correspond to the for- 
mation of embryos shown in the left, middle, and right 
columns of Figure [TJ From panel C it is evident that 
^^g®^ increases sharply each time the maximum gas tem- 
perature T^^^ exceeds the crystallization temperature 
Tcr ~ 800 K (horizontal dash-dotted line). After each 
such episode, prolonged periods of gradual decline in 
follow, which are caused by dilution of the disk material 
with pristine amorphous dust accreted onto the disk from 
the infalling envelope. 

In the last two episodes, T^^^ goes above the dust 
evaporation temperature Tevap and this leads to an im- 
mediate sharp drop in ^c.?- This exemplifies cases with 
''evap < '?"dest- The dust evaporation timesc ale can be 
approx imated (excluding factors of unity) as (|Navakshinl 
[2010i]) 

For the typical embryo masses of Memb=5-20 Mj, the 
corresponding Tevap lies in the 4000-60000 yr range. In 
fact, this is a lower estimate simply because embryos may 
experience several contraction and expansion episodes 
before actually reaching Tevap- Embryos that form in the 
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Fig. 1. — Gas surface density (g cm~^, top row), gas midplane temperature (Kelvin, middle row), and crystalline silicate fraction (bottom 
row) at three times since the formation of the central star: t = 0.074 Myr (left column), t = 0.132 Myr (middle column), and t = 0.178 Myr 
(right column). 



early disk evolution are of lower mass and closer to the 
star than those that form later in the embedded phase 
(because disk grows in mass and size with time). As 
a consequence, first embryos take longer time to evap- 
orate dust and shorter time to migrate into the star, 
an effect found in our nun i erical simulations and also 
noted by ICha fc NavakshinI (|2QlQf ) . For these first em- 
bryos, migration timescales are just several orbital peri- 
ods (Torb=700-1600 yr) and are often shorter than Tevap- 
As time passes and more massive embryos on wider or- 
bits start to form, the migration/destruction time scale 
may exceed that of dust evaporation, as indeed occurred 
during the last two episodes of crystalline dust produc- 
tion. 

Migration of embryos into the disk inner regions and 
through the sink cell, as manifested by strong 5 x 



10 ^ Mq yr ^) mass accretion bursts, also contributes 
to the decrease in (^^]f^ but leads to an increase in . 
Note that episodes 3 and 4 are not associated with strong 
mass accretion bursts, indicating that the corresponding 
embryos are dispersed by tidal torques exerted by spiral 
arms, before they actually fall through the sink cell onto 
the star. The process of disk fragmentation diminishes 
after t ~ 0.2 Myr, no embryos survive beyond this time, 
and the subsequent evolution shows a gradual decline in 
^^g®^. From this moment on, ^^.s^^ represents the crys- 
talline silicate fraction in the disk only. 

We do not take into account a possible formation of 
solid cores (and associated reduction in disk enrichment 
with crystalline dust) but we can estimate its efficiency 
based on the work of Navakshin (2010a). For Memb = 
3-10 Mj, Nayakshin found Tsed ^ (4-6) x 10^ yr. This is 



Crystalline silicate production 



Time (Myr) 

0.2 0.3 



Time (Myr) 

0.2 0.3 




o 



o 

OS 





-I— > 

OS 

o 

'U) 
(/) 

o 

Q. 

o 
E 
< 






Q. 

E 


-I—' 

C/) 

03 
D) 

E 
x 



2500 



2000 



1500 



1000 



maximum temperature 
crystalline silicate fraction 



500 




accretion rate onto the star 
crystalline silicate fraction 




7e-5 



6e-5 >s 

® 



5e-5 
4e-5 
3e-5 





c 
o 





o 

2e-5 o 

CO 
CO 



1e-5 



0.0 0.1 



0.2 



0.3 0.4 



0.0 0.1 0.2 0.3 0.4 0.5 



Time (Myr) 



Time (Myr) 



Fig. 2. — Crystalline (panel A) and amorphous (panel B) silicate fractions in the disk (black), envelope (green), and star+sink cell 
(blue) as a function of time passed since the formation of the central star. Panel C shows the maximum temperature in the disk (red) and 
crystalline silicate fraction in the disk (black), while panel D also shows the mass accretion rate onto the star (red). 

We have run another eight models with /3 = (3.2 — 12) x 
10"^ and Mcore = 0.3 - 1.25 M©. The net result is that 
the dust processing is correlated with disk fragmentation 
and models that do not form embryos (e.g., models with 
low Mcore and /3 and consequently with disks of low mass 
and size) show little crystalline silicates, reinforcing our 
conclusions. 

We stress that our proposed mechanism can coexist 
with other crystallization mechanisms, localized to the 
inner few AU, and complements them by providing a 
direct source of crystalHne silicates at large distances. 
Some outward transport may still takes place and ac- 
count for the fact that comet Wild 2 has crystalline sil- 
icates with the chemical and isotopic composition simi- 
lar to that of the chondri tic meteorites originate d in the 
terrestrial planet region (iMcKeega n et al.ll2QQ6f ). Con- 
versely, chondritic meteorites may have formed in the 
atmospheres of protoplanetary embryos via dust sedi- 
mentation and released into the inner few AU when the 
embryos were dispersed by tidal torques near the young 
Protosun. The apparent presence of refractory Calcium 
Aluminum Inclusions in the cometary material does not 
invalidate this scenario since temperatures in at least 
some of the embryos can exceed 1300 K (panel C in Fig- 
ure [2]). We do not take into account possible opacity 
dependence on the composition and properties of dust 
grains due to the complexity of physics involved. Nev- 
ertheless, dust growth and crystallization are likely to 



comparable to or less than Tevap 

> (4-16) X 10^ yr, which 
implies that solid cores may form before dust is vapor- 
ized. At the same time, r^ed is comparable to or longer 
than Tdest = a few x Torb ^ a few x 10^ yr, certainly 
for those embryos that form first in the early embed- 
ded phase and are characterized by shortest migration 
timescales. This simple timescale analysis suggests that 
the solid core formation is feasible and may reduce the 
efficiency of crystalline dust enrichment but is unlikely 
to shut down the proposed mechanism completely. The 
latest numerical hydrodynamics simulations seem to con- 
ffim that the solid core formation is not expected to work 
for every embryo in the disk (jCha fc Navakshinll2QlQ[ ). 

5. SUMMARY 

Thermal annealing of amorphous dust grains in the hot 
interiors of massive fragments (or planet/brown dwarf 
embryos) forming in protostellar disks during the early 
embedded stages of stellar evolution, followed by subse- 
quent dispersal of the fragments, is a promising gateway 
for the production of crystalline silicates. This mech- 
anism can release processed dust to various radial dis- 
tances from (sub-)AU scales, if fragments are destroyed 
on their approach to the central star, to hundred- AU 
scales, if they are dispersed by tidal torques exerted by 
spiral arms. More studies, including the effect of dust 
growth and solid core formation, are needed to further 
explore the efficiency of this mechanism. 
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decre ase opacities (|Mennella et al.|[T998l : iD'Alessio et al.l 
l2QQl[ ) , thus enhancing disk fragmentation and crystahine 
silicate production. 
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